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We study the role of the superconducting proximity effect on the electron-phonon energy exchange 
in diffusive normal metals (N) attached to superconductors (S). The proximity effect modifies the 
spectral response of the normal metal, in particular the local density of states. This leads to a 
weakening of the electron-phonon energy relaxation. We show that the effect is easily observable with 
modern thermometry methods, and predict that it can be tuned in structures connected to multiple 
superconductors by adjusting the phase difference between superconducting order parameters at the 
two NS interfaces. 



I. INTRODUCTION 

The state of an electron system subject to driving de- 
pends on three ingredients: on the external driving, on 
the response functions of the system, and on internal and 
external relaxation within the system and out of the sys- 
tem. For sufficiently strong but constant driving, the 
system can be brought out of equilibrium from its sur- 
roundings, in which case the steady state can be deter- 
mined from a balance between the driving and the relax- 
ation. By far, the most relevant relaxation mechanism 
for electrons in metals is caused by the coupling between 
electrons and phonons. If the dominant relaxation mech- 
anism is known to a high accuracy, the result of such heat 
balance can be used to study the driving. An example 
of such a procedure takes place in hot-electron thermal 
radiation detect or SfiiS, which rely on the changes in the 
electron temperature due to changes in the amount of 
radiation coupling to the device and driving the electron 
system. Typically such devices require accurate ther- 
mometry through an easily measurable observable that 
is sensitive to the electron temperature. One such fre- 
quently employed thermometer is based on a junction 
between normal metals and superconductors. 

When a normal metal is brought in contact with super- 
conductors, the superconducting order parameter leaks 
out to the normal side. This superconducting proximity 
effect"^ changes the spectral response of the system^ but 
it also affects the relaxation mechanisms. In this paper 
we study how this effect manifests itself in the electron- 
phonon coupling and consider the schematic structure 
shown in Fig. [TJ Some of the features induced by the 
proximity effect are the changes in the local density 
of states^ and a finite pair amplitude inside the nor- 
mal metal. For example, the density of states obtains 
a minigap whose size Eg{(j)) depends on the phase dif- 
ference 4> between the superconductors:^ at </> = for 
£^Th < A, « 3.12£;Th whereas Eg ^ aX 4> ^ t:. 
Here i?Th = TiD/L^ is the Thouless energy of the normal 
metal piece with a diffusion constant D and length L and 
A is the energy gap of the superconductor. 

By utilizing quasiclassical equations, we show how the 
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FIG. 1: Schematic structure considered in this paper: the 
superconducting electrodes (assumed to be large and resid- 
ing at the phonon temperature Tp) induce a proximity ef- 
fect into the normal metal. This proximity effect changes the 
spectrum of the electronic excitations and thereby also the 
electron-phonon energy exchange Qe-p- The spectrum can 
be controlled via the phase difference cf> of the superconduct- 
ing order parameter. This allows one to explore the effect as 
shown below in Fig. [6] 



collision integrals describing the electron-phonon cou- 
pling are changed by the proximity effect. We consider 
in particular the resulting effects on the electron-phonon 
relaxation rate, and the heat power flowing between the 
two systems. Our examples are calculated for a SNS sys- 
tem where the normal metal is sandwiched between two 
superconductors. In this case the proximity-induced ef- 
fects can also be tuned by adjusting the phase difference 
between the order parameters of the two superconducting 
contacts. 

This paper is organized as follows. At first we de- 
rive the general diffusive-limit electron-phonon collision 
integrals from Usadel equation^ for Keldysh Green's 
functions. Then we study especially the case when the 
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electron system can be described with Fermi function 
at a local temperature Te and the phonon system with 
Bose function at temperature Tp. For small deviations 
\Te — Tp\ <C Te the collisions can be described with a 
relaxation rate, whose dependence on energy, tempera- 
ture and (j) E^re shown for an example geometry. Then 
we concentrate on the energy current between the two 
systems, the quantity entering heat balance equations. 
This quantity is well defined for arbitrary temperature 
differences between the electron and phonon systems. Its 
dependence on position, length of the normal metal and 
(j) are exemplified for example geometries. Finally we 
suggest how the phase dependence can be measured in 
a system easily manufactured with present experimental 
techniques. 



II. THEORETICAL FRAMEWORK 

When describing nonequilibrium superconductivity, it 
is convenient to concentrate separately on two types of 
excitations^iS "longitudinal", described by the antisym- 
metric part /l = f[—E) — f{E) of the electron energy 
distribution function with respect to the chemical poten- 
tial of the superconductor, and "transverse", described 
by the symmetric part /t = 1 ~ f{E) — f{—E). The pre- 
vious describes "thermal" excitations whereas the latter 
is relevant for charged excitations. The detailed forms of 
the collision integrals for the two types of excitations are 
derived in the Appendix. At quasiequilibrium described 
by an effective temperature Tg, the most relevant type 
of excitations are the longitudinal ones. In this case, as- 
suming the elastic mean free path exceeds the phonon 
wavelengthfi the collision integral is given by 



J dujUjhgniuj)K{e,Lu)R{e,Lu;Te,Tp). 



(1) 



Here A is the electron-phonon coupling constant, the ker- 
nel 

K{e,Lu) = N{e)N{e + u;)^ F,ie)F,{e + cj) - Fsie)F,sie + uj) 

describes the changes in the spectrum of the junction, 
N{e) is the (local) reduced density of states, and Fc{e) 
and Fs{e) are projections of the pair amplitude (anoma- 
lous function) as described in Eq. ([9]). The part de- 
pending on the phonon and electron temperatures is in 
quasiequilibrium 
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The kernel K{e,u!) is in general position dependent, and 
thereby it makes also the collision integral depend on the 
position. Once this collision integral is known, it can be 
inserted in a kinetic equation, such as those presented in 
Ref. i. 



In what follows, we describe the electron-phonon scat- 
tering inside the SNS junction in terms of the scatter- 
ing rate and the heat current flowing between the elec- 
tron and phonon systems. In order to find the density of 
states and the anomalous function inside the normal re- 
gion of the SNS junction, we solve numerically the Usadel 
equatio: 



4.7.8.9 



(g«V3«) 



(2) 



where 7 is a small positive parameter describing inelastic 
scattering and the retarded Green's function satis- 
fies the normalization condition (g^)'^ = 1. We assume 
that the NS interfaces are clean and that g^ obtains the 

form of bulk Green's function g^^^iy_ = ej^/ J (\ — A^fs + 



A/y^e^ — A2[cos((/i)if2 +sin(0)jfi] at these boundaries. 

Here e+ = e -I- Z7, A is the absolute value of the super- 
conducting order parameter, </> is its phase and fi are the 
Pauli matrices in Nambu space. For lower-transparency 
junctions the proximity effect and thereby the effects de- 
scribed below will be reducedji 



III. SCATTERING RATE 

Close to equilibrium, the electron-phonon scattering 
can be described with an energy dependent scattering 
YSiie}^ The latter can be obtained by assuming that the 
Keldysh part of the scattering self-energy is related to the 
retarded (R) and advanced (A) parts via = {Yi^ — 
T,^)fL and similarly for the Green's function, g^ « (9^^ 
9^){fL + Sh). This gives us II « Tl^lSfL with 



Vl^l = -Tr[{g^-9^){^^'-^^)] 



A / dujijj^sgn{uj)K{e^Lj)- 
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sinh 
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(3) 



In the absence of superconductivity, N{e) « 1 is approxi- 
mately a constant and Fc^ Eg = 0. In this case at e = 
we get rf_p = 7AC(3)(fcsT)3^ 

The scattering rate is plotted in Fig. [5] as a function 
of energy for a few phases </> and at two different tem- 
peratures. The plotted rate is averaged across the weak 
link. 



and normalized to r^^_„ to illustrate the corrections due 
to the proximity effect. At energies e > fcsT the rate 
rises as in the absence of superconductivity. The main 
modification due to the proximity effect is the phase- 
dependent minigap, which causes a huge drop in the scat- 
tering rates. Moreover, above the minigap the rates are 
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slightly lower than in the absence of the proximity effect 
— this is due to the fact that within a thermal coher- 
ence length — v^i5/(27rfcj3r) from the NS interfaces 
electron-phonon coupling is weakened at any value of the 
phase (see Fig. [3]). 
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FIG. 2: (color online) Phase and energy dependence of 
the electron-phonon scattering rate Fe-p normalized to the 
normal-state zero-energy result F^_p = 7AC(3)(A;Br)^ at two 
different temperatures (top: ksT — iSxh/S and bottom: 
ksT — S-Etii)- We have assumed Ao = 20iJTh and aver- 
aged the rates across the normal-metal piece. The solid line 
shows the energy dependence obtained in the normal case for 
which K{e,uj) = 1. Note the logarithmic scale of the figure. 
The exact low-energy values in the presence of the proximity 
effect depend on the magnitude of the inelastic scattering pa- 
rameter 7 in Eq. ([2]). Here and throughout we have chosen 
7 = IQ-^Et^. 



IV. HEAT CURRENT 

The most relevant quantity in describing the electron- 
phonon interaction is the heat current associated with a 
temperature difference in the electron and phonon sys- 
tems. This quantity can be described also far from equi- 
librium, i.e., for arbitrary difference between the electron 
and phonon temperatures. The heat current density is 



obtained by multiplying 1^ with energy e and the normal- 
state density of states vp at Fermi energy, and integrating 
over the energy: 



P.-r 



Xvp j d€dujeuj'^sgn{uj)K{e,u})R{e,uj;Te,Tp). 



(4) 



In the normal state, the integral can be evaluated ana- 
lytically, with the resuli^i^ 



(5) 



where E — 24C(5)i^i?fc^A. Measured values of E in differ- 
ent metals are tabulated in Ref. fl!. 

In most cases, the heat diffusion within the normal 
metal is much stronger than out from it. In this case the 
variations in the temperature of the normal metal become 
small and the average temperature is determined by the 
rate integrated over the volume of the junction, i.e.. 



(6) 



Here A is the cross section of the normal metal. In the 
normal case is independent of position and thus 

Q^_p — P^pfl, where V, = AL is the volume of the 
normal metal. 

Proximity effect changes the electron-phonon heat cur- 
rent drastically. Most of the changes are limited w ithin 
a thermal coherence length — \/T^D/ {2TTkBT) from 
the superconductors. This is illustrated in Fig. [3l which 
shows Pe-p normalized to P^_p as a function of posi- 
tion and temperature. Note that changing Tg inside the 
normal metal does not affect the superconducting gap as 
the superconductors are assumed to be at the phonon 
temperature. 

A relevant factor for the SNS system is the relation of 
the energy gap inside the superconductor to the Thouless 
energy of the weak link. This ratio equals to the square 
of the ratio between the length of the junction and the 
superconducting coherence length. Figure H] shows the 
position-averaged heat current as a function of electron 
temperature for different lengths of the junction. At low 
temperatures below the minigap, ksT ^ Eg, the heat 
current is vanishingly small because of the lack of states 
inside the normal metal. However, also at higher tem- 
peratures there is some reduction, whose measurement is 
well within reach of state-of-the-art thermometry. 

Apart from the temperature and length of the junc- 
tion, the nature of the proximity effect is quite sensitive 
to the phase difference (j) between the superconductors. 
This sensitivity is mainly due to the phase-dependent 
minigap in the density of states. This is illustrated in 
Fig. O where Qe-p is plotted as a function of (j)- The 
strongest dependence takes place naturally at tempera- 
tures close to the phase-dependent minigap. However, 
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FIG. 3: (color online) Position and temperature dependent 
heat current density Pc-p between the electron and phonon 
systems inside the normal part of the SNS system of length L. 
Here = 0, A = 20-ETh, Tp = and the heat current density 
has been normalized to the normal-state value = ST^ . 








FIG. 4: (color online) Dependence of the heat current on elec- 
tron temperature and the length L of the junction compared 
to the zero-temperature coherence length ~ ^^hD/A. Here 
(?!> = and Tp = 0. 

some phase dependence is still remaining even at higher 
temperatures. 

In the long-junction regime L ^ = and temper- 
atures up to some IQEx/kB the temperature dependence 
of the power Qc-p can be fitted rather well to a curve of 
the form 

Qe-p = Qe%CXp(-T7T), (7) 

where T* — cE'Yii/kB is a temperature scale describing 
the proximity-induced variations. For the prefactor we 
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FIG. 5: (color online) Phase-dependent electron-phonon heat 
current at a few different temperatures with Tp = and A = 
20i?Th. The closing minigap for — > tt shows up as a rapid 
increase of Qo-p at low temperatures. For fcsTe somewhat 
larger than the minigap, the phase dependence becomes much 
weaker although the overall heat current is still smaller than 
in the normal state. 

get c ~ 3.5 . . . 3.7, depending on how large temperatures 
are included in the fit. This behavior resembles that of 
bulk superconductorSfi^ but where the gap A is replaced 
by T* . However, in the SNS case the exponential behav- 
ior extends to comparatively larger temperatures than in 
bulk superconductors where one has to assume fc^T ^ A 
for the law to be valid. 

As seen in Fig. [SJ Qe-ph becomes independent of the 
phase at temperatures larger than the zero-phase mini- 
gap. In this case the temperature dependence can again 
be expressed with Eq. However, for low temperatures 
the effective temperature scale becomes of the order of 
the phase-dependent minigap Eg{(j)). 

V. EXPERIMENTAL DETERMINATION 

A possible experimental setup to measure the effects 
discussed above is shown in Fig. [HJa) . The structure can 
be fabricated through standard lithography techniqueSfi 
The device consists of a radio frequency superconducting 
quantum interference device (SQUID) ^lii where a su- 
perconducting loop is interrupted by a N wire of length 
L. The strength of the electron-phonon interaction in 
N can be modulated periodically by an externally ap- 
plied magnetic field, which gives rise to a total flux <i> 
through the loop area. Neglecting the inductance of the 
superconducting loop, the phase difference becomes then 
(f) — ^/ (/i/2e). The SQUID allows magneto-electric char- 
acterization of the SNS junction, and thus also the de- 
termination of some of the relevant parameters of the N 
wire. For instance, the Thouless energy of the junction 
can be extracted from the temperature dependence of the 
SQUID critical supercurrent. Alternatively, the phase 
can be controlled by an externally driven current in the 
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case when the superconductors are not connected. How- 
ever, in this cases only phases (j> « — 7r/2. . .7r/2 can be 
accessed. The N region is connected to four additional su- 
perconducting electrodes through oxide barriers, so to re- 
alize normal metal-insulator-superconductor (NIS) tun- 
nel junctions. The NIS junctions are used to heat (or 
eventually to cool) the electrons in N, and as sensitive 
thermometers to measure T^.^ Both the NS contacts and 
the NIS junctions provide nearly ideal thermal isolation 
of the N region and therefore in the following we neglect 
the thermal conductance through the superconductors. 
One further advantage of using tunnel junctions directly 
connected to N stems from the fact that the state of the 
electron system and the superconducting correlations in- 
duced by proximity from S electrodes will be virtually un- 
affected by the presence of tunnel-coupled probes.^* The 
basic requirements for this are that the currents driven 
through the NIS junctions arc much smaller than the SNS 
critical current, and that the tunnel junction resistances 
are much larger than the resistance of the normal metal. 
These requirements are easy to satisfy in practice. 

Figure [DJb) shows a sketch of the relevant thermal 
model of the wire in the proposed setup. Upon heat- 
ing electrons with a constant power Qh provided by the 
NIS heaters the steady-state Te established in N will de- 
pend on the energy relaxation mechanisms occurring in 
the wire. In metals at low lattice temperature (typically 
below 1 K), the main relaxation mechanism is related to 
electron-phonon interaction which in the present setup 
is strongly phase dependent. In practice, what will be 
measured is the power flowing between the electron and 
phonon systems averaged over the wire volume fl, i.e., 
Oc-pCTc, Tp, 0). At fixed lattice temperature the pre- 
dicted steady-state phase-dependent thus follows from 
the solution of the thermal-balance equation 

QH + Qc-p(Te,Tp,0) = O. (8) 

The resulting electron temperature Te can then be 
probed with the NIS thermometers. 

The result of a such a heat balance is plotted in Fig. [5] 
which displays the steady-state electron temperature Tg 
for some values of Qh as the phase is varied and the bath 
temperature is Tp = 50 mK. To illustrate the strength of 
the effect, we chose a typical E^^/kB — 100 mK. For 
each Qh the electronic temperature is strongly modu- 
lated by the phase, and it decreases by an increasing (p 
due to the enhancement of the electron-phonon interac- 
tion for phase difference close to tt (see Fig. 5). In par- 
ticular, for Qh — 0.1Sr2(i<^Th/fcB)^, variation of Te of 
the order of 40 mK (0.4 E'xh) can be achieved by vary- 
ing the phase (flux), while even for the large injected 
power Qh = 100Sr2(i<^Th/fcs)^ the variation amplitude 
is somewhat below 30 mK (0.3 £^Th)- 

Let us discuss some of the typical materials parameters 
for the proposed measurement. By choosing, for instance, 
aluminum (Al) as S material (with Aq = 200 /ieV) and 
copper (Cu) as N wire (with D = 0.01 m^s~^ and E = 
2 X 10^ Wm-3K-5, see Ref. T) and length L = 0.8 /xm, 
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FIG. 6: (color online) (a) Scheme of a possible experimental 
setup. A radio-frequency superconducting quantum interfer- 
ence device (rf-SQUID) containing a SNS junction is threaded 
by a magnetic flux "1> which allows to tune the electron-phonon 
interaction in the N region. Additional superconducting elec- 
trodes tunnel-coupled to the N wire allow to probe the aver- 
aged phase-dependent Qc-p{Te,Tp,(j>). These serve both as 
heaters (h) and thermometers (th). (b) Sketch of the ther- 
mal model of the N wire (see text), (c) Steady state elec- 
tron temperature inside the SNS junction as a function of 
the phase for a few driving powers at a phonon temperature 
Tp = £JTh/(2fcs). For this curve, we chose iJTh/fcs — 100 mK 
(see text). 



we get Ao/Et = 19 and Erh/kB = 10 fieV = 120 
mK. Moreover, if we assume that the wire is 50 nm 
thick and 400 nm wide, the "unit power" in Fig. [5] is 
T.n{ETh/kB)^ = 0.8 fW. 



VI. CONCLUSIONS 

In bulk superconductors, the presence of the energy 
gap leads to a strong suppression of the heat transport 
between the electron and phonon systems. In this pa- 
per we show that a similar effect can be expected for nor- 
mal metals in close proximity to the superconductors and 
suggest how it can be probed. Additionally, the heat cur- 
rent between the two systems can be controlled in situ by 
changing the phase difference between two superconduct- 
ing contacts. The suppression of the heat conductance is 
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especially relevant in thermal radiation detectors — for 
example in the Proximity Josephson Sensor (PJS) sug- 
gested in Ref.H. This is because of two factors: First, the 
weaker the thermal relaxation, the stronger is the tem- 
perature rise for a given power of radiation. Secondly, 
the main cause of noise in these devices is related to the 
strongest thermal relaxation channel, which typically is 
the electron-phonon coupling. Reducing this coupling 
will therefore also reduce the noise. Both of these factors 
will improve the sensitivity of such thermal detectors. 



Here A — — is the spectral function and g^'^'^ 
are the retarded, advanced and Keldysh Green's func- 
tions, respectively. From these self-energies, the collision 
integral appearing in Usadel equation is 
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We are mostly interested in the Keldysh part, which 



has two components: II 
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/2. In the following we parameterize the 



Green's function by g — g h — hg , h = JlTq + frTs, 
and 



VII. APPENDIX 



9 



-9 



The self-energy for the electron-phonon scattering in 
the case of thermal phonons at temperature Tp and with 
the general Keldysh Green's function for the electrons is 
given byi^ 



where 

E«/-4 ^ A 



5^c— p — 
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-uj'^sgn{uj) 



—LU sgn(w) 



E^ 



coth 



UJ 



2kBTp 
i(e + cj) 



coth 



iksTp 



Moreover, the normalization condition implies g^ — — 
1. The Advanced Green's function can then be obtained 
by g^ = —'T'3{g^)^T3- Note that inside a proximity struc- 
ture, all the three variables, g, f and (fi are in general 
complex-valued. Moreover, we define N — Ke[g] (local 
density of states) and the real quantities 



==Re[/cos((/.)], Fs 
=Im[/cos(,^)], F, 



Re[/sin(0)] 
Im[/ sin{(j))]. 
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After some lengthy but straightforward algebra, the col- 
lision integrals become 



It =A / dLuujhgniuj){ [Fc{e + uj)Fs{e) - F,(e + uj)Fc{e)] 
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+ [N{e)N{e + uj)+Fcie + Lu)F^{e)+Fsie + u)Fsie)]coth 
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and 



/i=A / dLoujhgTi{u)\ [N{e)N{e + u)-F,{e)F,{e + uj)-Fsie)Fsie + u)]x 



coth 



/ LU 



\ 2kBTp 
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+ [Fs{e + Lo)Fc{e) - Fc{e + uj)Fs{e)] 
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+ [Nie)Nie + u;) ~ F,ie)F,ie + lj) - F,(e)F,(e + uj)]hie)hie + to) 



These simphfy considerably in the quasiequihbrium hmit where /t = and = tanh[e/(2fcBT)]. In that case one 
obtains the coUision integral presented in Eq. Ilj. 
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